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1. Introduction 

The zero-range process is a stochastic lattice gas where the particles hop randomly with an 
on-site interaction and the jump rates g(n) depend only on the local particle number n. It 
was introduced in [HI as a mathematical model for interacting diffusing particles, and since 
then has been applied in a large variety of contexts, often under different names, (see e.g. 
and references therein). The model is simple enough for the steady state to factorize, on the 
other hand it exhibits an interesting condensation transition under certain conditions. When 
the particle density exceeds a critical value p c , the system phase separates into a homogeneous 
background with density p c (the fluid phase) and all the excess mass concentrates on a single 
lattice site (the condensate). 

Besides spatial inhomogeneities (see e.g. flU), condensation can be caused by an effective 
attraction between the particles if the jump rates g{n) have a decreasing tail as n — > oo. A 
generic class of such models with a power law decay 

g{n) ~ 1 + b/n 1 as n — > oo , (1) 

with positive interaction parameters 6,7 has been introduced in [0 O, and condensation 
occurs if < 7 < 1 and b > 0, or if 7 = 1 and b > 2. Results on homogeneous zero- 
range condensation have been applied to many clustering phenomena in complex systems 
such as network rewiring Q, traffic flow (HI or shaken granular media (9l[T0l, for a review 
see 0. Using a mapping to exclusion models, zero-range condensation can also be used as 
a generic criterion for phase separation in driven diffusive systems with one or more particle 
species [11]. The condensation transition in this model is now well understood lTT2l [T3l . 
also on a mathematically rigorous level lfT4l [T51 [T6ll . and many variants have been studied 
El \T]\ \M EE! US ETJ- In [231 the influence of specific non-random perturbations have 
been studied for models with asymptotically vanishing jump rates. 

The assumption of strict spatial homogeneity is not very realistic in applications to real 
complex systems which often exhibit disorder due to local imperfections. In ll24ll a randomly 
perturbed version of the model CD has been introduced, and it turned out that an arbitrary 
small perturbation has a drastic effect on the critical behaviour. Using heuristic arguments 
it was shown that condensation occurs only if < 7 < \, significantly changing the phase 
diagram of the unperturbed system. These first results only applied on finite systems and 
crucial questions on the distribution of the critical density and whether or not the system 
exhibits condensation in the thermodynamic limit remained open. In this paper, we provide 
rigorous results on the quenched free energy and prove the existence of a finite critical density 
in the thermodynamic limit. We also give accurate expansion results to compute the value of 
thermodynamical variables and their distributions for small and large perturbations which are 
supported by detailed simulation results. 

The paper is organized as follows: In Section [2] we define the perturbed zero-range 
process and introduce thermodynamic variables of interest (such as the free energy) and our 
numerical methods. In Section [3] we derive rigorous results in the thermodynamic limit and 
provide expansion results for small and large noise in Sections [4] and [5] In Section [6] we 
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conclude and discuss the relevance of the thermodynamic limit results in real finite systems 
using exact numerics for fundamental diagrams from recursion relations. 

2. Definitions and numerical methods 

2.1. The disordered zero-range process 

We consider a regular, <i-dimensional lattice A of finite size |A| = L with periodic boundary 
conditions. A configuration is denoted by (r) x ) x <=A where rj x E {0, 1, . . .} is the occupation 
number at site x. The dynamics of the zero-range process is defined in continuous time, such 
that with rate g x (f]x) site x E A loses a particle, which moves to a randomly chosen target site 
y according to some translation invariant probability distribution p(y — x) . For example in one 
dimension with nearest neighbour hopping, the particle moves to the right with probability p 
and to the left with I— p. Our results do not depend on the specific choice of p(y — x) as long 
as it is of finite range. For simplicity of presentation, we focus on jump rates g x given by 

9x (n) = e ^ {n)+b/n '' for n > 1 , g(0) = . (2) 

Here a, b and 7 are positive parameters and x E A, n E N are independent, identically 

distributed random variables with 

E[£ x (n)] = 0, E[^.(n) 2 ] = l, and 5 := logE[ e -^ (ri) ] < 00 . (3) 

By Jensen's inequality and strict concavity of the logarithm we have 5 > 0. For a = 
the asymptotic behaviour of the jump rates is given by (OQ) so the present model, which has 
been introduced in ll24l . can be interpreted as a perturbation of the generic homogeneous 
model. Note that the exponential form of the jump rates and the use of standardized random 
variables £ x (n) is purely for notational convenience, since jump rates have to be non-negative. 
In simulations, we will mostly use uniform random variables £ x (n) as a generic case, but 
also discuss the simplifications in the case of Gaussians in Section 14. 1[ Our results hold 
for any sufficiently regular, generic perturbation of © which does not change the expected 
asymptotic behaviour of the rates, i.e. for which Hindoo E[g z (n)] = g < 00 is finite and the 
same for each lattice site x. Perturbations with x-dependent asymptotic expectations lead to 
spatial inhomogeneities which act as an additional, independent source of condensation. This 
mechanism is different from the one induced by the asymptotic decay of the jump rates which 
we want to focus on in this paper, and has been studied previously in simpler models (see e.g. 

mi 

The main difference to these studies on purely spatial disorder [4] is that in the present 
model the noise also depends on the occupation number at each site. In [|24l it was shown that 
this feature leads to a drastic change in the phase diagram for finite systems with L fixed. As 
opposed to the unperturbed model CQ), the perturbed model © exhibits condensation only for 
7 E (0, 1/2). This is due to the contribution of the perturbation to the partition function, and is 
formulated precisely in Proposition Q] in Section [3l A particularly suitable application where 
such perturbations are relevant is for example the bus route model [12311 . where (QQ) describes 
the rates at which a bus proceeds to the next stop given distance n to the previous bus. The 
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larger the distance, the slower the rate due to more passangers queueing at the stations, which 
can lead to a condensation of buses (on a route with periodic boundary conditions). To test 
robustness of this phenomenon in realistic situations, it is very natural to assume that the 
actual functional behaviour of the rates is of the form © with a small random dependence on 
the bus x. 



2.2. Stationary distributions and thermodynamic variables 

It is well known (see e.g. Il26l 13) that the above zero-range process has grand-canonical 
factorized steady states v 1 ^ = ® x eA^- The single-site marginals are 

^( n ) = — 7 — c w x {n) for n > , (4) 
%x (A* J 

where the chemical potential fixes the particle density, and the stationary weights w x 

are given by the jump rates via 



w x {n) = 

k=l 



n 

\\gxik)- 1 = exp ( - aS x (n) - /3(n)) . (5) 



Here the contribution from the perturbation, S x (n) = J2k=i^(k), can be interpreted as 
the position of a random walk on R after n steps with independent increments £ x (k). The 
contribution from the interaction is f3(n) = X^=i^A 7 ' an d acts as an n-dependent drift 
term. This holds independently of the translation invariant jump distribution p(y — x) and for 
each realization of the £ x (n), i.e. iA is a quenched distribution. The single-site normalization 
is given by the partition function 

oo 

z *(/-0 = y^ ex P f n A* - vSxjn) - /3{n)j (6) 

n=0 

which is strictly increasing and convex in fx. For convergence of z x (fj) < oo it is necessary 
that the drift term (3{n) dominates the stochastic part oS x (n). This has been used in Il24ll to 
predict the phase diagram of the model, which we review in Section [3] together with our new 
results. 

We denote 

f x (ji) =logz x (/J l ) , (7) 
and the local density can be calculated as the derivative 

Px[p) ■= {Vx)v M = d , (8) 

and it is a strictly increasing function of ji. Here (•}„ denotes an average over the quenched 
distribution u^, and z x , p x are still random variables w.r.t. the perturbation. Averages w.r.t. 
realizations of the £ x (ri) are denoted by E and determine the behaviour of the system in the 
thermodynamic limit. As usual for disordered systems, we consider two quantities. The 
annealed free energy is defined as 

fM:=logE[z x (fi)}, (9) 
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and it can be interpreted as the free energy of a site with an average perturbation. Using © it 
can be rewritten as a shift of the homogeneous free energy 

oo 

f A (p) = log ex P ( n V + n5 ~ P( n )) = f(P + 6 ) 7 (10) 

n=0 

where /(//) and z(p) (and in general all quantities without subscript) refer to the homogeneous 
system with a — 0. Given that z(p) is monotone increasing in p, /(//) is as well, and hence 
/a (a*) > /(/•*) f° r ai l Z 1 an d a U o" > 0. On the other hand, the quenched free energy 

f Q (p):=E[logz x (p)]=E[f x (p)], (11) 

is the average of f x and is the physically relevant quantity in the thermodynamic limit. 
In general, by Jensen's inequality and concavity of the logarithm the annealed free energy 
/a (a*) > /q(aO provides an upper bound. 

As is clear from the tail behaviour of © (see Section [3] for a precise statement), the 
grand-canonical partition function z x (p) does not exist for p > 0, and evaluating the above 
quantities at p = yields the critical values of thermodynamic system parameters [24]. In 
particular, the (average) critical density is 

p c :=E[p x (0)} e [0,oo] , (12) 

which we prove to be finite if and only if 7 G (0, 1/2) in Theorem [3j This is the main rigorous 
result of this paper. 



2.3. Numerics 

Using (flOl) . the annealed free energy can be computed to arbitrary precision from the 
unperturbed systems. The issue is to generate reliable numerical results for quenched 
quantities under the presence of disorder. Convergence of z x (p) (and hence of f x {p)) is 
assured by the analytical results shown in Section[3l however, at the critical point the effect of 
the perturbation is maximal and convergence can be very slow. Therefore it is useful to define 
and analyze truncated quantities 

N 
n=0 

and, analogously, f£ and . By definition, 

f x {p) f x {p) a.s. (almost surely) , 

i.e. convergence holds for all typical realizations. But the speed of convergence is random 
and depends on x. As usual, we will use empirical averages on finite systems with |A| = L 

El[/^(m)]:=^E^) (14) 

to estimate E[f^(p)]. Provided that the latter is finite, the law of large numbers implies 
convergence of E^f^} as L — > 00 for all N > 0. Furthermore, since f x and E[/^] are both 
monotone increasing in N, we have 

f Q (p) = lim E L \f»fr)] , 
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Figure 1. Asymptotic behaviour of \f£ (0)] for 6=1 and several values of a with uniform 
distribution of £, x (n). (left) Fixed L = 10 4 with increasing N plotted against 1/N. As 
expected, E L [/^(0)] diverges for 7 = 0.6 > 1/2 (dashed lines), and for 7 = 0.4 < 1/2 
(full lines) it increases monotonically to an upper bound which is the numerical estimate of 
/q(0). Note that for7 < l/2the smaller a, the faster the convergence, (right) Fixed N = 10 5 
with increasing L for 7 = 0.4, the error bars represent the standard error of the mean. For 
small L, rare large values dominate the empirical average at the critical point, and have to be 
compensated choosing L large enough. 



and the limits for N and L commute. 

We find that for generic values b « 1 and 7 < 0.4 we get reliable estimates of for 
all jj, < using 

N = 10 5 and L = 10 4 , (15) 

provided perturbation strengths are a < 1. This is illustrated in Fig. Q] for fj, — 0, where it is 
shown that approximations become largely independent of the truncation parameter N and the 
number of samples L above those values. We use these parameters and a uniform distribution 
of the £, x (ri) in all simulations in the paper, except explicitly stated otherwise. Analogous 
estimates hold for p c and E[p x (//)] with a < 0.8, since convergence turns out to be slower for 
this observable (not shown). 

In general, convergence becomes slower as a increases. A rough estimate for the minimal 
truncation parameter N with fi — is the point where the drift contribution /3(n) in © 
exceeds the random walk contribution S x (n) and thus successive terms of the sum start 
decreasing. For large n, 

/3(n) ~ — - — n 1 ' 1 and \aSJn)\ ~ a^/n , 
1-7 

which leads to N ~ a 1 - 2 '/ (cf. also Section^. This can grow very fast with a, in particular for 
7 close to 1 /2, and numerical results for large perturbations are computationally expensive. 
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3. Rigourous results in the thermodynamic limit 

3.1. Preliminary results 

Before we state our main new result we summarize some simple facts for completeness, some 
of which have already been used in Il24l . 

Proposition 1 Let a > 0. z x (p) and p x (p) are a.s. smooth functions of p for p < 0. For 

7 G [1/2,1] and allb > we have 

z x (/i) < oo a.s. p < and z x (p) — > oo a.s. as p /* . (16) 

For 7 G (0, 1/2) and allb > we have 

z x (p) < oo a.s. <^ p < and Zx(/i) — > z x (0) < oo a.s. as p /* .(17) 

77?e same statements hold for p x (p) and all higher moments, in particular, for the site 
dependent critical density 

p x {0)<ooa.s. <=> 7 G (0,1/2). (18) 
Proof. The law of the iterated logarithm (see e.g. lETTl ) implies for the random walk part of 

limsup } aSx ( n ) J= = aV2 a.s.. (19) 
n^oo vnloglogn 

Furthermore, we have 

n , , 

Mn) = V — ~ n l ~"t as n ->• oo , (20) 

v ' ^ kri 1 - 7 
fc=i ' 

and thus /3(n) ^> juS^n)! a.s. as n — >■ oo if and only if 7 G (0, 1/2). This implies 

- aS x (n) - c/3(n) ->■ -00 a.s. and n i e -°s*{n)-cfi(n) a.s. (21) 

for all c > 0, g > 0. In particular, w x (n) — > a.s. as n — > 00, and convergence is fast 
enough to bound the sum z x (0) = X] n >o <EH> implies that 

P(w a; (?7,) > e _/3( - n ^ 2 for only finitely many n) = 1 , 

and therefore there exists a C > such that 

Zx(0) <CJ2 e~^ n)/2 < 00 a.s. . (22) 

n>0 

The same holds for higher moments, e.g. for the density we have 

F(n w x (n) > e _/3( - n ) for only finitely many n) = 1 , 

which, together with (1221) . implies that p x (0) < 00 a.s. . If 7 G [1/2, 1) this argument only 
holds for u < 0. For p = and 7 G [1/2, 1) (fT9l) and (1201 imply that the random walk part 
dominates w x (n) for large n. Therefore 

F(w x (n) > 1 for infinitely many n) = 1 
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On fixed finite systems with lattice A this directly implies that for 7 6 (0,1/2) 

pM) := xJ2 p*(°) < 00 a - s - ' (23) 

i.e. there exists a finite critical density, which depends on the realization of the perturbation. 
In correspondence to previous results such as [15] the process is then expected to exhibit 
condensation when the number of particles diverges, which can be formulated in terms of 
canonical distributions and the equivalence of ensembles (cf. Section [6]). In the following we 
focus on properties of the grand canonical distributions in the thermodynamic limit. 

Proposition 2 The annealed free energy ([9]) is /a(/0 = /(/^ + S) E [0, 00], given in terms of 
the unperturbed model for all /iGi For all 7 £ (0, 1) and for 7 = 1 and b > 1 we have 

f A (p) < 00 & H<-&. (24) 

Proof. By direct computation we have 

00 00 

E[zM\ ^EK(n)] =^E[e- ff «-«] n e-^ + ' in = zfa + S) (25) 

n=0 n=0 

by monotone convergence, since the terms in the sum are positive. The rest follows immedi- 
ately from the well-known properties of /(//) = log z(fi) for the unperturbed model (see e.g. 

EU2!). □ 

This implies in particular, that even though z x (0) < 00 a.s. for 7 £ (0, 1/2) and it is 
a well defined random variable, we have E[^.(0)] = 00. However, log z x (0), which is the 
thermodynamically relevant quantity, has finite expectation /q(0) in that case, as we will see 
in the next subsection. 

3.2. Main results 

Theorem 3 For the quenched free energy 4771) we have 

f(ji)<fQ(ji)<f(ti + S) = f A (ji) farallpeR. (26) 
For 7 £ (0, 1/2) we have 

f Q (ji) < 00 fi < and p c = E[p x (0)] < 00 . (27) 

For 7 £ [1/2, 1] we have 

/g(/i) < 00 n < and /g(/i) — > 00 as fi /*■ , (28) 

anJ p c = 00. 
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Figure 2. Upper and lower bounds d26l i for the quenched free energy Jq (full line) given 
by the annealed and unperturbed free energies Ja (dashed lines) and / (dash-dotted lines), 
respectively. Parameters are b = 1, a = 0.3 (implying S = 0.0446) for uniformly distributed 
noise £ x (n), and 7 = 0.4 and 0.6. /q diverges for 7 = 0.6 only very close to fi = (cf. also 
Fig. [U, which cannot be inferred from the behaviour of upper or lower bounds. Data for /q 
are simulated according to ( fT4b and fTBV 



By the definition of the critical density for finite systems (1231) . and the fact that the summands 
are independent positive random variables, we have a strong law of large numbers 

p c (A) — > p c a.s. as L — > 00 , (29) 

which covers condensation with p c < 00 for 7 G (0, 1/2), as well as the case p c = 00 for 
7 G [1/2, 1]. Fig. [2] illustrates the bounds (|26|) on the quenched free energy, which are rather 
accurate for small perturbation strength a away from the critical point, but do not contain any 
information on whether the system condenses or not. 



Proof of Theorem 3. The upper bound in d26l) follows immediately from Jensen's inequality 
and Prop. [2l The lower bound follows from (1371) in Section HI 

Let 7 G (0, 1/2). In the following we prove /q(/u) < 00 for /i = 0, which implies the 
same for all ji < by monotonicity. Write z^(0) = J2n=o w x( n )- Then z®(0) = 1 and 



<(0) > 1 a.s.. Let A N 



ELCT) - f?-\0)\. Then 



EL/f +1 (o) 



/"(0)]=E[log(^ +1 (0)/4 V (0))] 



E 



log 1 



w x (N + l) 
(0) 



< E 



since log(l + u) < u for u > 0. With w x (N + 1} 
independence of ^ x (iV + 1) we get 



E 



w x (N 



(0) 



e «- 6 /(jv +1 r E 



w,(iV) 
^(0) 



(30) 



w :c (^)e-^( 7V+1 )" f '/( JV+1 ) 7 and by 



(3D 
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Now, w x (N) < z^ (0) a.s. and we can estimate 



E 



Since z£ (0) > 1 a.s. we have, using the asymptotic form of f3(n) given in (1201) . 

p (5w > 1/iy2 ) - PK(A ° > 1/n2) - ¥ {~ aSxiN) > \i^ Nl ~ 1 - 2l0 % N 

for all N large enough. Since 7 G (0, 1/2) this is bounded above by exp [— Af 1_27 (y3^) 2 (g^)], 
using standard estimates from moderate large deviations (see e.g. OT1 [321 and 031 for a 
general reference). Thus with (|3Q|) and (I3TT) we have Xlfcli ^fc < 00 an d 

TV 

/q(0) = lim E[/^(0)] = lim ^ A fc < 00 . (33) 

TV— >oo TV— s-oo iI — ' 

fe=l 

The second statement, p c < 00, can be shown very similarly. Let 

1 * 

(0) := -ttt nw *( n ) < 00 ( 34 ) 

x ^ ' n=0 

for all N > 0, which is well defined since 1/^(0) G (0,1] a.s.. With Ajv = p%{0) 
we get 

7V^(iV+l) 



< A" a.s. . 



E[Ajv] < iVp(^M > 1/iV 3 ) + 1/N 2 



^(0) 

Analogous to (|32|) this leads to 

' w x (N ) 

and since z x (0) > 1 a.s. we have 

Kttst > 1/iy3 ) - p(Wx(Ar) > 1/iy3) - F (- aS ^ > \x^ Nl ~" ~ 3logN ) 

for all N large enough. This is again bounded by exp [ — A rl ~ 27 (jz^) 2 (g^)] since the leading 
order term is unchanged. The rest follows analogously. 

For 7 G [1/2, 1], on the other hand, (1281) follows immediately from the almost sure 
behaviour of z x (p) given in (TT6l) . and for « < one can use simple exponential bounds anal- 
ogous to the above. □ 



This result implies that for 7 G (0, 1/2) the local particle density p x (0) at the critical 
point has finite mean. The same can be shown analogously also for all higher moments of 
the occupation number, which are given by higher order derivatives of f x (0). We will see in 
the next section using a heuristic expansion, that for small noise f x can be approximated as 
a small perturbation of / for the unperturbed system. The actual distribution of f x (0) is very 
hard to describe analytically or access numerically with adequate precision, since it has heavy 
(sub-exponential) tails as implied by the following result. 
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Proposition 4 Let 7 G (0,1/2) and b,a > 0. Then we have for all X > 

E[ e A ^(°)] = 00 , (35) 
i.e. the distribution of f x (0) does not have exponential moments. 

Before we proceed with the proof, we introduce some notation which is used again later in 
Section [47TI We would like to stress the dependence on a of the partition function and write 

00 

z x (p,a) = z x (p) = ^>/-<^W-/%) . 

n=0 

Recalling that z(fi) = Y^=o e n ^~^ n ^ for the unperturbed system, we can define a probability 
distribution p(n), n > where the random variable X takes the value — S x (n) with 
probability p(n) = g ~^y~ • Denoting (-) p as the expectation w.r.t. this distribution, we can 
write 

z x (p,a) =z(n)(e (TX ) p . (36) 

Proof of Proposition SI Writing f x (fi, a) analogously to z x (/i, a), we have 
E [ e A/*(o,<7)] = E[z x (0,a) x ] > E[z x (0,a)] X = 00 

for A > 1, directly by Jensen's inequality and (T25T) . For A G (0, 1) we use (1361) . and writing z x 
formally as an expectation we can apply Jensen's inequality to get 

z x (0, a) x = z(0) x (e° x ) x p > z(0) x (e x ° x ) p = ^O)^ 1 z x (0, Xa) . 
Now, taking expectation w.r.t. the disroder 

E[z x (0,a) x ] > z(0) x ^E[z x (0,Xa)} = 00 , 
which follows from (1251) and using that z(0) G (0, 00) for the unperturbed system. □ 



4. Expansion for small perturbation 

4.1. Free energy distribution 

The effects of the noise on the critical density and other thermodynamic variables are hard to 
quantify in general beyond the results in Section [3j but expanding the partition function for 
small noise a — > leads to reasonable approximations in comparison with the unperturbed 
system. With Prop. [T] we have for 7 < 1/2, ji < and all a > 0, 

z x (p, a) = z(fi) (e aX ) p < 00 a.s. , 

where we use the notation introduced in (1361) . Taking logarithms on both sides and using 
Jensen's inequality w.r.t. the distribution p we get 

log z x (p, a) = log z(/i) +log(e x ) p > hgz(p) + (X) p . 
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Figure 3. Convergence of the rescaled probability density function of f x for vanishing a, 
with \x = 0, b = 1, 7 = 0.4 and TV = 10 5 . For each value of er (symbols as in legend) 10 6 
samples were taken, (left) For uniform £ x in) the limit is non-Gaussian, and can be determined 
numerically from (l39i l as shown by the full line, (right) For Gaussian £, x (n) the limit is also 
Gaussian and can be determined explicitly (1401) . The data for larger a is skewed towards 
positive values due to higher order corrections. 



Under the expectation w.r.t. the realizations of S x (n) the second term vanishes and we obtain 
a lower bound for Jq, 

/qO*)>/0*)+E[Wp] = /(a*)- (3V) 
Expanding the exponential in (1361) around a = yields for a fixed realization of the noise 

z x { l M,a) = z{ l i)Y,— ] (*"V (38) 
z — ' ml 

m=0 

According to Prop. [T] all moments of z x (fi,a) are finite, so (X m ) p < oo a.s. for all m. 
Furthermore, \(X m ) p \ < C J2^Lo Pnn m ^ 2 is almost surely bounded by the m/2-th moment 
of p n , which is bounded by ml for all p < 0. Therefore the series (1381) converges absolutely, 
and z x (fi, a) (as well as f x (n, a)) are in fact analytic functions in a around a = for all fixed 
// < 0. 

From (1381) we get 

oo n 

log fefri, a)) = log (*(//)) + log (l + ^( Xn )p) > 
and expanding the second logarithm on the right hand side up to second order in a yields 

2 

f x (fi, a) = /(//) + a(X) p + °- ((X 2 ) p - (X)fy + 0(a 3 ) . 
Using S x (n) = Ylk=i £x(&) we get for the first order term 

oo oo oo 

Wp = - E ^( n >w = - E^ n ) -••• = - E &(*o*U*o > 

n=0 n=l fc=l 



Condensation in randomly perturbed zero-range processes 



13 



where F^{k) = ^2 n>k p{n) 1S ^ e ta ^ °f m e distribution p(n). Therefore, the deviation of 
f x (p, u) from the unperturbed system's free energy f(/S) scales with a, and 

-+-ituk)FM ■ (39) 

This is a sum of i.i.d. random variables with vanishing prefactors, which depends in general 
on the distribution of the £ x (k) (see Fig. [3]). In the particular case that the £ x (k) are Gaussian 
this is a linear combination of independent Gaussians, and therefore, 

f x (fj.,a) - f(p) 



N(0, 1) with s(p) = J2 F t*( k ) 2 < 00 > ( 4 °) 

fe>i 

and the fluctuations are Gaussian to leading order. 
4.2. Expected values 

Taking expectations in the expansion (1381) one can also estimate /q = E[f x ] and a) = 
E[p x (p, a)]. We have E[5 , a; (ra)] = 0, and also E[S x (ni)S x (n 2 )S x (n 3 )] = and so on for all 
odd powers. Therefore only even powers contribute and the expansion of the quenched free 
energy is 

2 

f Q (p,a) = E[logz x (^a)) = /(/i) + y (E [(X 2 ) p ]-E [(X) 2 p ] )+0(a 4 ) .(41) 

After some straightforward calculations summarized in the appendix the expectations can 
again be expressed in terms of the tails F^ik) = J2 n >kP( n )- This leads to expressions for the 
quenched free energy and the density 

f Q (fi,a) =/( At ) + a 2 0(/i) + O(a 4 ) 

E\p x fa a)] = p(fM) + a^'{p) + 0(a A ) , (42) 



where the coefficients are given by 

- oo 



2 

n=0 

Note that unlike the expressions for the distribution of f x (p), these results do not depend on 
the distribution of the £ x (n). As is shown in Fig. [4] the expansion coincides very well with 
numerical data for values of a even up to 0.5. 

For larger values of a higher orders contribute to the expansion. In a first attempt to 
understand these contributions one can compute the fourth order term for a 4 of the cumulant 
expansion, which is given by 



24 



-6(X)l+12(X)l(X\-3(XX-MX) p (X\+(X% 



Condensation in randomly perturbed zero-range processes 14 




Free energy Critical density 



Figure 4. Convergence of the expansions of the free energy /q(0, <r) (left) and the critical 
density p c = E[p x (0, a)] (right) as given in ( f42l > and ( |43T > for small disorder. Rescaled data for 
different values of 7 collapse well onto the predicted a 2 behaviour (full line) for values of a up 
to approximately 0.5. For a > 1 clearly higher order terms become relevant. Here N = 10 5 , 
L = 10 5 and errors are of the size of the symbols. 



In order to evaluate terms of the form ~E[S x (ni) S x {n 2 ) S x (n 3 ) S x (n4)} we can use nested 
conditional expectations and obtain 

E[S x (n 1 )S x (n 2 )S x (n 3 )S x (n 4 )} = 3n\ + m(n 2 - n x ) + ni(n 3 - n 2 ) = 2n\ + n x n 3 . 

It is difficult to find a simple formula similar to (1431) for the fourth order coefficient, but it can 
be computed numerically to arbitrary precision. However, including this term does not give 
any substantial improvement of the prediction for large a (not shown). Apparently, the higher 
order coefficients do not decay fast enough, and in particular for a > 1 the behaviour cannot 
be understood by looking only at the first few terms in the expansion. A different approach to 
tackle large disorder is presented in the next section. 



5. Large disorder 

For large a, the disorder term can dominate the exponent in w x (n) © up to relatively large 
values of n. For typical realizations with non- vanishing probability we can have 

'a(l - 7 )v / 2\ 2 /( 1 - 2 ^ 



crV2^ - /3(n) > as long as n < n* := h ' j , (44) 

where we have used the asymptotic form of (3{n) ~ jz^n 1 ~" / . The dominant contributions 
to the free energy in this case come from typical disorder realizations with large aS x (n) ~ 
a\/2n. Note that — > 00 with a — > 00, and it will serve as a large parameter in the following. 
Approximating the sum z x (0) © as an integral we get with u = n/n* replacing aS x (n) by 



POO 

z x (0) ^n* e nl ~~' c{u) du where 
Jo 
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Free energy Critical density 



Figure 5. Leading order behaviour of the free energy /q(0, a) (left) and the critical density 
p c (right) for large disorder. The numerical values (symbols as in legend) are expected to 
converge to the predictions (full lines) for a — > oo, and the leading order power in a is already 
well predicted also for the computationally accessible values for a < 10. Here N = 10 5 , 
L = 10 4 and errors are of the size of the symbols. 



c (u) =n'l~ 1 / 2 aV2u —u 1 ^ . (45) 

1-7 

The integral is dominated by the saddle point value u s where the exponent takes its maximum, 
i.e. 

c '( M ) = nr 1/2 4=- 6M ~ 7 = ' ( 46 ) 

\/2u 



so that 



. an:~ 1/2 \ 2 /(i-27) / 1 \ 2/(1-27) 
u s = t=— = 777; r < 1 • (47) 



V2b J V2(l- 7 ) 
Expanding c(u) m c(u s ) + \c"{u s ){u — u s ) 2 around the saddle point we get 



E [.o g ^(0)] « log (n, l—^e^) = - log^j + C( b ,;) 



1 — 



as ex 00. Here C(b, 7) = 2 _1 - 2 t b 1 - 2 ^ (1 — 27) /(l — 7) and the lower order terms grow 
only logarithmically in o and can be ignored for large a. 
A similar argument works for the critical density 



Pc 



00 00 



(48) 



n=0 n=0 

where both sums are dominated by the saddle point u s , leading to 

/ a \ 2/(1-27) 

p c w n*u s = (^-y^J • (49) 

Again, there are similar additive corrections as for log 2^(0) which can be ignored for large 
a, and we do not write them explicitly. Both predictions are confirmed relatively well by 
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numerical data (see Fig. 7). Since we can only go up to a = 10 due to numerical restrictions, 
there are still relatively large finite size corrections. But the asymptotic slope of the curves in 
a double-logarithmic plot corresponding to the leading order powers in a are well confirmed. 
Note that the corrections are smaller for small 7, since here the width of the Gaussian in the 
saddle point approximation proportional to l/(n;J;~ 7 |c"(u s )|) is smaller and the integrand is 
concentrated more sharply around the saddle point. 

6. Discussion 

In this paper we have provided a fairly complete picture of the influence of a generic 
perturbation on the condensation transition in zero-range processes with decreasing jump 
rates. Our results include a rigorous analysis of the grand-canonical measures and the 
associated thermodynamic quantities such as the free energy and the critical density. We 
also provide detailed numerical data to illustrate our results, and heuristic arguments to 
approximately predict the behaviour of the system for small and large disorder. 

In order to understand the relevance of our results for condensation in real systems we 
consider the canonical stationary distribution 7Ta,jv of a system with a fixed number of N 
particles on the lattice A = {1, . . . , L} with periodic boundary conditions. 7Ta,jv can be written 
as a conditional distribution 7Ta,at = ( . | Y^xeA^x = N), which is actually independent of 
p (see e.g. Ifl4l0 . For simplicity, we focus the discussion on totally asymmetric jumps with 
p(z) = 5 Zt i. In this case the average stationary current 

3a,n ■= (9(Vx)) WAN = Z ^' N ~ 1 , (50) 

is given by a ratio of canonical partition functions which can be computed exactly via the 
recursion relation (cf. e.g. ED EE!) 

N 

Za,N = 2^ W x( k ) Z A\{x},N-k • (51) 
k=0 

As a reminder, for the unperturbed model (OQ) we have z(p) < oo if and only if p < and 
p c < oo as long as 7 G (0, 1) or 7 = 1 and b > 2. In the thermodynamic limlit N, L — > 00, 
N/L — > p, it has been shown (see e.g. [51 0) that the canonical distributions converge to a 
grand canonical factorized distribution with spatially homogeneous marginals as given in 
©. For p < p c , p is chosen to fix the density p(p) = p via relation ([8]), and for p > p c , the 
maximal possible value p = (corresponding to density p c ) is chosen independently of p. In 
the latter case the system phase separates into a homogeneous background at density p c with 
distribution vq, and a condensate where a macroscopic fraction of all particles concentrates 
on a single lattice site [5][T6]|. A particularly useful signature of the condensation transition 
is the behaviour of the expected current {g{rj x )) v in the thermodynamic limit as a function 
of the density. It is monotonically increasing for densities below p c and becomes constant 
for densities above. Convergence to this limit is typically accompanied by particularly strong 
finite size effects which have been studied in OTTl . The current shows a characteristic non- 




condensing fluid 

Figure 6. Canonical current-density relation (150b for a condensing system with 7 = 0.4 (left) 
and a non-condensing system with 7 = 0.7 (right). For small disorder both systems appear to 
be condensed, and the thermodynamic limit result is only valid for large disorder or extremely 
large system sizes which are not accessible numerically. 

monotonic behaviour, consisting of an increasing fluid branch and a decreasing condensed 
branch, as can be seen in Fig. [6] for a = 0. 

For the perturbed model we have shown in this paper that the parameter range for which 
condensation occurs in the thermodynamic limit changes to 7 G (0, 1/2). Nevertheless, for 
small noise a and 7 below as well as above the critical value 1/2 the current shows the same 
characteristic behaviour, and the system appears to be condensing as shown in Fig. |6j Only 
for rather large values of a > 1 the system for 7 > 1/2 (Fig. [6l right) appears fluid for all 
densities, as the analysis of the grand-canonical measures predicts. The large fluctuations in 
the current for intermediate densities result from changes in the condensate location due to 
the environment. The behaviour shown in Fig. [6] for system sizes L = 500 is typical for all 
numerically accessible sizes up to L = 5000. The thermodynamic limit results do therefore 
not give a good approximation of the behaviour of finite systems with moderately large system 
sizes, which are particularly important in many applications, such as shaken granular media 
ll9l[T0l or traffic flow M- 

While the most interesting properties of the site-dependent free energy f x , such as finite 
mean and sub-exponential tails, are included in our results, it would be interesting to estimate 
the exact tail behaviour of its distribution. This requires an understanding of the leading 
order contributions to the partition function which is an interesting question in itself. While 
properties of exponential functionals of Brownian motions with constant drift are known 
to great detail (see e.g. OOlO . the form of the weights w x (n) in the present model do not 
allow for an exact analysis. First numerical results indicate a crossover in the behaviour, 
where depending on the system parameters the sum is dominated by a large number of small 
contributions or a small number of large contributions. 
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Appendix. Calculation of expansion coefficients 

In the following we compute the coefficients (1431) of the expansion in Section 14.21 Since S x 
has independent increments, ~E[S x (n)S x (k)] = min{n, k}. The first expected value in the 
bracket in (HTT) is Y^=o n P( n ) = Pit 1 )- F° r me second expected value we need to compute 
terms of the form 

oo oo oo 

min{n, k] p{n)p{k) = — np(n) 2 + 2 np(n)F^(n) , 

n,k=0 n=0 n=0 

where F^(k) = J2 n >kP( n )- Now usep(n) = F^n) — F^(n + 1) and the trick 

(F» - F,(n + = \{F,{nf - F,(n + l) 2 + (F» - F> + l)) 2 ) . 

p(n) 2 

Then summation by parts 

oo oo 

^Tn(F> + l) 2 -F>) 2 ) = -]TF,( 

n=0 n=0 

leads to 



n? 



min{n, k} p(n)p(k) = F f 

n,k=0 n=0 

One finally obtains 



oo oo 



e [(x\ - (x)l] = p (jm, o)-J2 F M 2 = E - ^W) • 



n=0 n=0 



An estimate of the density follows by differentiation of the free energy expansion w.r.t. /i. It 
is useful to compute first 



<9u 

n>/c n>fc 
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